On the low-frequency spatial dispersion in wire media 
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The work is dedicated to the theoretic analysis of wire media, i.e. lattices of perfectly conducting 
wires comprised of two or three doubly periodic arrays of parallel wires which are orthogonal to 
one another. An analytical method based on local field approach is used. The explicit dispersion 
equations are presented and studied. A possibility to introduce a dielectric permittivity is discussed. 
The theory is validated by comparison with the numerical data available in the literature. 
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I. INTRODUCTION 

In the recent years the periodic metallic lattices have 
found many applications both in optical and microwave 
ranges (see, for example, in 1 and 2 ). However, some fun- 
damental problems have not been resolved yet, even for 
typical metallic electromagnetic crystals. One of them is 
the problem of low-frequency spatial dispersion in wire 
media (WM). The low- frequency spatial dispersion of a 
simple wire medium (a doubly periodic regular array of 
parallel wires) has been studied only recently in 3 . In the 
present paper this theory is generalized for double and 
triple wire media. The study of spatial dispersion effects 
in the above mentioned variants of WM has been started 
in work 4 . However, this study (based on the numerical 
approach) is far from being complete. Our theory sig- 
nificantly complements the results 4 . It is analytical one, 
and in order to validate it, a comparison to the results 
from 4 is carried out. 

The unit cells of lattices under study are shown in Fig. 
1. They are comprised of two (2d or double wire medium) 
or three (3d or triple wire medium) doubly periodic reg- 
ular arrays of parallel infinite wires which are orthogonal 
to one another. The wires are assumed to be perfectly 
conducting. The host medium is a uniform lossless di- 
electric with permittivity £o and permeability fi . We 
denote the radii of wires directed along x,y and z-axes 
as r x , r y and r z , respectively. The periods of the lattice 
along x,y and z-axes are denoted as a, b and c, respec- 
tively. The lattices are spatially shifted with respect to 
each other by half period (see Fig. 1). The wires axes 
positions in the chosen coordinate system are determined 
by equations: 

• the x-directed wires: y = bn + b/2 and z = cl + c/2, 

• the y-directed wires: x — am 4- a/2 and z = cl, 

• the z-directed wires: x — am and y = bn, 

where m, n and / are integers. 

In order to model an electromagnetic response of a wire 
we apply the local field approach. We assume that the 
wires diameters are small as compared to the wavelength. 
Thus, every wire can be described in terms of effective 




FIG. 1: Unit cells of double wire medium (a) and triple wire 
medium (b). 



linear current referred to the wire axis. The wire with 
radius r oriented along a unit vector d (|d| = 1) can be 
characterized by a " polarizability" a, relating the com- 
plex amplitude I of the induced current and the local 
electric field E loc : 



7 = a(r ,A:,q-d)E loc -d. 



(1) 



Here k is the wave number of the host medium, and qd = 
gii is the longitudinal component of the wave vector of 
propagating mode. The following expression for a was 
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obtained in 5 : 



a(r ,k,q\\) 
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(2) 



where 7 w 0.5772 is the Euler constant and 77 = /^o/eo 
is the wave impedance of the host medium. 

Let the eigenmode under consideration have the wave 
vector q = (q x ,q y ,q z ) T . Expressing the local field pro- 
duced by all wires except the reference one through the 
current induced in the reference wire we obtain the dis- 
persion equation. It relates the components of q with k 
(i.e. with frequency u>). Then we can introduce effec- 
tive material parameters of the wire medium which fit 
this dispersion equation. In this paper we widely use the 
results obtained in our precedent papers 3 ' 5 for a simple 
WM. Therefore, in the next section a very short overview 
of those works is presented. 



II. SIMPLE WIRE MEDIA 

The geometry of a simple wire medium comprising the 
z— directed wires is shown in Fig. 2. 
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FIG. 2: Simple wire media: a doubly periodic lattice of par- 
allel ideally conducting thin wires. 

The currents in the wire with numbers (ra, n) (counted 
along x and y axes, respectively) are related with current 
/ induced in the reference (zeroth) wire through the wave 
vector q: 



I„ 



Jgjilxam+qybn) 



(3) 



Since the field re-radiated by the wire (m, n) is propor- 
tional to I m . n , the local electric field acting to the zeroth 
wire can be expressed in terms of the so-called dynamic 
interaction constant C: 



E l ° c = C(k, q x , q y , q z ,a, b)I, 
where (see in 3 ' 5 ): 



(4) 



C(k 1 q x ,q y ,q z ,a,b) 



V (k 2 -q 2 z ) 
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(5) 



Rm,n = \J (am) 2 + (bn) 2 and all (m, n) except m — n = 
are summed up. The expression (5) can be rewritten 
in the following form 5 : 



C(k,q x ,q y ,q z ,a,b) 



E 



V(k 2 ~ ql) 
2jkb 

■ 1 (n) 
sin/Ox a 



1 



• ; (0) 

sin k x a 



1 (0) ; (0) 

Kx cos k x a — cos q x a 



~7q \kx cosk x n ^a — cosq x a 27r|n|^ 
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Where 




(0) 



(7) 



and we choose Re{y / (J} > 0. Those formulae physi- 
cally correspond to the representation of the WM as a 
set of parallel grids (of z— directed wires) located paral- 
lel to one another with period a along the a;— axis (see 
also Fig. 4 for the case of 2d WM). Every grid radiates 
the spectrum of Floquet harmonics with wave vectors 
(k x , q y + 2irn/b, q z ). The series with summation over n 
in the right side of (6) describes the contribution of the 
high-order Floquet modes into electromagnetic interac- 
tion of those grids. The dispersion equation follows from 
(1) and (4): 

[a _1 (r 0) A;,g z ) - C(k,q x ,q y ,q z ,a,b)] 1 = (8) 

Taking into account expressions (2) and (6) one can 
rewrite (8) in the following form: 



(k 2 - ql) 



1 1 b 



1 



sin k 



(°), 



bk x cos k x a — cos q x a 



1 



• ; (n) 
sin k x a 



1 



\bk x n ^ cosk x n ^a — cos q x a 27r|n| J 
This equation has three types of solutions: 



1 = 0. (9) 
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1. Ordinary waves, in case when I = in (9). They 
have no electric field component along wires (E z = 
0) and propagate without interaction with the lat- 
tice. Their dispersion plot corresponds to the host 
medium and is shown in Fig. 3 by thin lines. 

2. Extraordinary waves, in case when the expression 
in square brackets in (9) equals to zero. They cor- 
respond to the nonzero currents 1^0 and have 
the nonzero longitudinal component of electric field 
E z =^ 0. Their dispersion properties are described 
in details in 5 ; their dispersion curves are presented 
in Fig. 3 by thick lines. 

3. Transmission-Line Modes (TLM), in case when 
(k 2 — q 2 ) = in (9). Those waves propagate along 
the wires, they are TEM waves (E z = 0), but / ^ 0. 
Their dispersion equation q 2 = k 2 has no restric- 
tion for components q x ,q y , and the phase shift of 
the currents in the adjacent wires can be arbitrary 3 . 




r x 
Normalized wave vector, qa/n 
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FIG. 3: Dispersion curves of wire media with filling ratio 
/ = Tir^/a 2 — 0.001 (square lattice). Thin lines - ordinary 
waves, thick lines - extraordinary waves. TLM are not shown. 

Under the quasi-static limit ka <C 2-7T and \q\a <C 27T, 
the dispersion equation for extraordinary waves trans- 
forms to 



r = ix + <$ + il = k ' 2 - k n, 

where the following notations are used: 

2ir/s 2 



log 



27T7-0 



+ F(r) 



(10) 



(11) 



s = vab, r = a/b and F(r) = F{l/r) is given by 

, 1 , ^ / coth(7rw) - 1\ 
F(r) = -^ogr + J2{ '-^ J + 



irr 
~6' 



(12) 



Parameter fco corresponds to the effective plasma fre- 
quency of the lattice loq = fc /\/£oMo- ^ or sc l uare lattices 
a = b one has F(l) — 0.5275. Comparing (10) with the 
well-known dispersion equation of uniaxial dielectrics we 
obtain an effective relative permittivity e of Id WM in 
the following form: 



e = ez z + x x + y yo, 



e(k,q z ) = 1 



k 2 



k 2 



(13) 



(14) 



The dependence of dielectric permittivity on q z given by 
(14) does not disappear untill frequency becomes zero. 
That fact means the low-frequency spatial dispersion of 
wire media. There is no low-frequency spatial dispersion 
for the extraordinary waves in the only case when the 
wave propagates across the wires (q z =0). At low fre- 
quencies the propagation of those waves can be described 
in terms of plasma-like permittivity e = 1 — k 2 /k 2 (sec 
also 6 ). Relat ive to those waves the wire medium behaves 
as a cold non-magnetized plasma (a continuous dielec- 
tric medium). In other propagation directions the wire 
medium behaves differently. In 3 we discuss the impor- 
tance of the low-frequency spatial dispersion in Id wire 
media. Below we theoretically show this phenomenon in 
2d and 3d WM. 



III. DOUBLE WIRE MEDIA 

Now, let us consider a double wire medium which is 
comprised by y— directed and z— directed wires. It is 
shown in Fig. 4 as a set of parallel grids located along the 
x— axis. Below the local field approach is going to be ap- 
plied taking the same approximation as it was done in our 
work 7 (where we have studied a doubly negative meta- 
material in a similar way). Thus, the approximation is 
following: the electromagnetic field produced by a single 
grid of wires at the distance from the grid a/2 is consid- 
ered as a field of a sheet of the average current J. That 
approximation is accurate enough under the condition 
when the wavelength in the matrix is large as compared 
to the grid periods (kb <C 27r and kc <C 27r) and period a 
is not smaller than periods b, c. In this case the y-oriented 
grids interact with z-oriented grids by the fundamental 
Floquet harmonic. Other harmonics are evanescent and 
their contribution into this cross-polarized interaction is 
negligible. 

We can express the (m, n)-numbered z— directed cur- 
rent through the reference (zeroth) z— directed current 

The same rule holds for y— directed currents: 

4 m ' l) (z) = I ye i{i^ m +i^+ivv) . (16) 
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FIG. 4: The structure of a double wire medium is represented 
as a set of planar wire grids. 



It should be noticed that the expressions in brackets in 
(23) are exactly the dispersion equations for the simple 
wire media (from y- wires and z— wires, respectively). 

In (23) coefficients C yViZZ are defined from (6), (19), 
(20). Now, let us calculate coefficients C zy and C yz us- 
ing the approximation of current sheets which has been 
mentioned above. The z— component of the electric field 
produced by a sheet with surface current J y (y,z) = 
yo(Iy/c)e^ qvV+q ^ z ' > at the arbitrary distance x from the 
grid can be expressed by the formula 8 : 

E v<z (x,y,z) = r mij Z!y {y,z)e-^\, (24) 



k x = -j\lql + q 2 z - k 2 . 



Summing up (24) over all to— numbered layers we can 
write: 



+ OC 



The currents I z and I y are related to the local electric 
fields acting on the z— and y— directed reference wires 
E l ° c and E y ° c through polarizabilities a z ,y' 

I z = a z E l ° c , I y = a y E y oc , a y , z = a(r y ^ z ,k,q y ^ z ), 

(17) 

which fit (2). Both E l ° c and E y oc contain contributions 
of z— and y— arrays: 



k 1oc = fS z ~> -+- fS v ^ f, 1oc = fSv) 4- fM 



(18) 



Co-polarized terms E z z ^ and E y could be expressed 
through I z and I y applying (4): 

E^=C ZZ I Z , C zz = C{k,q x ,q y ,q z ,a,b), (19) 



c zy =m- > 



jq^am p -j\m-\ |fe x a 



2 kx kc 

The summation result is following 



(25) 



,-, _ VQyQz je lq * a/2 cos(q x a/2)sm(k x a/2) 

u zy - ~m — ~ i, : ■ \ M ) 



k<-r kc 



cos q x a — cos k x a 



Taking into account the phase shift between z and y grids 
we obtain 

C yz — £~ lqw C zy c/b 



Wh je 3q ^ a/2 cos (ggfl/2) sin (k x a/2) 
k x kb cos q x a — cos k x a 



(27) 



E i = C vyIy> C yy = C(k,q x ,~q z ,q y ,a,c). (20) 

(v) (z) 

The cross-components E z and E y could be ex- 
pressed through I y and I z , respectively: 

E^=C Z yIy, E^=CyJ Z . (21) 

The cross-polarized interaction factors C yz ^ zy are evalu- 
ated below. Substituting equations (19), (20) and (21) 
into (17) we obtain a system of equations 



{Cyy - a y 1 )7„ + C yz I z = 
C zy I y + {C zz - a" 1 )!, = 



(22) 



First of all, it should be noticed that the solution of 
(22) when I ViZ = corresponds to the ordinary waves 
with polarization along the a;— axis, propagating in the 
plane (y — z). The dispersion equation for such waves is 
q 2 + q 2 — k 2 , q x — 0. Setting the determinant of (22) 
equal to zero we obtain a dispersion equation for the ex- 
traordinary waves (I y , z ^ 0): 

{Cyy - oTy X ){C zz - O - C yz C zy = 0. (23) 



Substituting (2), (6), (26) and (27) into (23) we derive 
an explicit dispersion equation: 



E 



1 1 c 

-log 



1 sin k. r a 



ir 2nr y ck x cos k x a — cos q x a 



smp y 'a 



\cf3 y n ^ cos/3y™^a — cos q x a 27r|n|^ 



(fc 2 - «2) 



1 1 b 

-log 



1 sinfc^a 



7T 2iTr z bk x cos k x a — cos q x a 



1 sm/3 z a 



tf V bp z n) cos f3 z n) a - cos q x a 27r|n| 



M^qI ( cos (q x a/2) sin (k x a/2) 



(28) 



k 2 bc \ cos q x a — cos k x a 



■5 



where 



IV. TRIPLE WIRE MEDIA 




The signs of all square roots are chosen so that 
R C {VT)} > 0- The dispersion equation (28) cannot be 
simplified (except the quasi-static limit) even in a spe- 
cial case when r y — r z and a = b = c. In fact, the per- 
fect square can be obtained in the left side of (28) if one 
neglects the contribution of high-order Floquct harmon- 
ics expressed by n— series. However, this approximation 
leads to the wrong results for the shape of isofrcquencies. 
Therefore we do not use it. 

The preliminary analysis of (28) reveals some special 
solutions. There are two solutions which correspond to 
TLM: the first is q y = k, q z = 0, q x is arbitrary; the 
second one is q z = k, q y = 0, q x is arbitrary. Those 
waves propagate either along the y— wires (when the elec- 
tric field averaged over the lattice unit cell is polarized 
along z) or along the z— wires (when the averaged elec- 
tric field is polarized along y). They are TEM waves as 
well as TLM in simple WM 3 . The component q x is a 
free parameter for TLM and plays a role of a phase shift 
between the currents in the adjacent grids of wires 3 . At 
the first sight, it seems strange that the electric field with 
non-zero z— component can propagate along y across the 
z— directed wires below the "plasma" frequency (which 
is the cut-off frequency for such waves in Id WM). How- 
ever, it is possible. When the TLM propagates along the 
y— wires, all grids of z— wires are excited, however the 
superposition of their fields exactly vanish in the planes 
x = am + a/2, where the grids of y— wires are located. 
This result can be easily obtained analytically for arbi- 
trary non-zero q x and q t = k. Due to the same reason it 
is also possible for y— polarized TLM to propagate along 
z. 

When q x = ir/a or q y q z = the right side of (28) 
equals to zero and the equation splits into two separate 
equations similar to (9) and describing the extraordinary 
waves in two simple WM. For q y = (or q z = 0) the 
absence of the interaction between two simple WM is 
trivial since the propagation holds in the plane (x — z) 
(or (x — y)) and the electric field is polarized orthogonally 
to y— directed wires (or to z— wires). However, the inter- 
action between two Id WM is also absent when q x = it /a. 
At low frequencies ka < 1 the equation q x — it /a corre- 
sponds to the excitation of TLM in both y- and z-arrays 
with polarization directions alternating along x. The ex- 
istence of this kind of TLM (which does not transport 
energy at all) is specific for 2d WM. 

More detailed study of (28) requires numerical calcu- 
lations and their results are presented below. 



Analysis of triple wire media can be carried out in the 
same way as it was described above. Similarly to (22) we 
obtain: 



(C xx a x )I y + C xy I y + C XZ I Z — 

C yX I X + (Cyy — Cty )I y + C yZ I Z = 

C ZX I X + C zy I y + {C zz - a, 1 )!, = 



(29) 



where C 
(6) and 



yy,zz,yz are determined by (19), (20), (26), (27) and 



C x 



C(k 1 q y ,q z ,q Xl d yi d z ). 



(30) 



Here on is denoted as a, = a^r,, fc, <?») and the subscript 
i means the Cartesian components (x,y, z). Other inter- 
action factors are as follows: 



C xy j 



nq x q y cos (q z c/2) sin (k z c/2) ^ 
k z ka cos q z c — cos k z c 



jqz c/2 



c = . yqxQz cos (q y b/2) sin (k y b/2) ^ qyb/2 
xz k y ka cos q y b — cos k y b 



c _ . VQxqy cos (g z c/2) sin (fc z c/2) e _ jqzC / 2 
yx k z kb cos q z c — cos k z c 



C z 



, Tjg x g z cos (q y b/2) sin (k y b/2) e _ jqyb/2 
k y kc cos q y b — cos k y b 



(31) 
(32) 

(33) 
(34) 



-J 



qi 



fc 2 , 



i\hl + iv - fc2 - 



There are no ordinary waves in that medium since there 
are no vectors orthogonal to all wires simultaneously. De- 
terminant of (29) gives the dispersion equation: 

(C xx - a- l ){C vy a^){C zz a-') (35) 



-{C xx — a x 1 )C yz C zy (C, 



mi 



a y X )C XZ C Z 



(Czz Q z *)C xy C yx -\- C xy C yz C zx + C xz C zy C yx 



0. 



The dispersion equation (35) with substitutions 
(19),(20),(30)-(34) is the final result for the triple wire 
media. The explicit equation is cumbersome and cannot 
be simplified in the general propagation case. However, 
in the special case when q x = all cross-polarized in- 
teraction terms (31)-(34) vanish, and the system (29) 
splits into two separate sets: the first one is the disper- 
sion equation of the Id WM C xx — a^ 1 , the second one 
is the system (22). The first case corresponds to the ex- 
traordinary waves propagating normally to the x-wires 
without interaction with y— and z— wires. There is no 
spatial dispersion for those waves (see above). The sec- 
ond case corresponds to the in-plane propagation in 2d 
WM, which will be studied below. In the present paper 
we do not consider the general case of the wave propaga- 
tion in 3d WM. 
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V. DISPERSION DIAGRAMS AND 
ISOFREQUENCIES OF A DOUBLE WIRE 
MEDIUM 

The dispersion diagram of a double WM for the in- 
plane propagation (q x = 0) of the extraordinary waves 
obtained by numerical solution of (28) is shown in Fig. 5. 
The chosen parameters of the wire lattice are a = b = c, 
r y = r z . The filling ratio is / = 2nrl/a 2 = 0.002. We 



1.6 




Normalized wave vector, qa/n 



FIG. 5: Dispersion diagram of a double wire media with filling 
ratio / = — 0.002 (cubical cell and equal radii). Thin 

lines - modes of the host medium (singular points of equation 
(28)), thick lines - modes of the 2d WM. 

use notations T = (0,0, 0) T , Z = (0,0,tt/c) t , and L = 
(0, n/b, ir/c) for the central point, the z— bound point 
and the corner point of the fundamental Brillouin zone, 
respectively. 

One can notice the significant difference between Fig. 
5 and the dispersion diagram of a simple wire medium 
(see Fig. 3). In the Fig. 5 one can see within the inter- 
val L — r two extraordinary modes which do not vanish 
at low frequencies k < ko and are not TLM. In simple 
WM the waves with nonzero longitudinal (with respect to 
the wires) component of the electric field cannot propa- 
gate at low frequencies since the phase shifts between the 
adjacent wires are small and the re-radiation of parallel 
wires suppresses the wave. In 2d WM it becomes pos- 
sible due to the electromagnetic interaction of the two 
orthogonal wire arrays. This is the result of the cross- 
polarized interaction of wire arrays. There are terms in 
C yz that cancel out the terms in C yy and C zz which are 
responsible for the suppression of the waves propagating 
obliquely in simple WM at low frequencies. 

The horizontal lines ka/(2Tt) = 0.1 and ka/(2n) = 0.3 
in Fig. 5 correspond to isofrcquency contours presented 
in Fig. 6 and 7, respectively. The isofrequency con- 
tour located around the L point is very unusual (close 
to the hyperbolic one). In Fig. 6 one can see that the 
contours of isofrequencies are rather close to four asymp- 
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FIG. 6: Isofrequency contours for double wire media at 
ka/(2-w) = 0.1. Two cases q x = and q x = 7r/(2a) are pre- 
sented. 



totes q VtZ = ±k. In spite of the rather low frequency as 
compared to u>o, the isofrequency contour located around 
the r point (q = 0) basically differs from the isofrequency 
of an isotropic dielectric (a circle) . Only in a special case 
of the in-plane propagation the isofrequency centered at 
the T point has practically circular shape and the phase 
velocity of this mode coincides with that of the host 
medium. When q x ^ the shape of this isofrequency 
becomes super-quadric and modes with hyperbolic isofre- 
quency tend to the same asymptotes q y _ z = ±k. When 
Qx = 7r/a the isofrequencies coincide with the asymptotes 
exactly. This case corresponds to TLM discussed above 
(which do not transport energy). The plot in Fig. 6 in- 
dicates the possibility of the two refracted waves (both 
extraordinary waves) for the rather large sheer of inci- 
dence angles. This effect keeps at the quasi-static limit. 

The 2d WM with proper orientation of wires with 
respect to the medium interface can possess the low- 
frequency negative refraction. It follows from the fact 
that the angles between the group and the phase veloc- 
ities for the mode corresponding to the hyperbolic con- 
tours in Fig. 6 and Fig. 7 can be close to 7r/2 (the 
normal to the isofrequency contour shows the direction 
of the group velocity vector) . 

At the frequencies close to the plasma frequency wo 
and higher two other modes appear with isofrequencies 
centered at T. They are shaped as two crossing ellipses. 
The modes with isofrequency curves close to q y>z = ±k 
are still present. The isofrequency contours for such a 
case (corresponding to ka/2i: = 0.3, q x = and q x = 
n/(2a)) are shown in Fig. 7. When q x increases at fixed 
frequency, the hyperbolic isofrequency contours in the 
plane (q y — q z ) approach the asymptotes in the same way 
as it happens for lower frequencies. The elliptic contours 
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FIG. 7: Isofrequency contours for double wire media at 
ka/(2n) = 0.3. Two cases q x = and q x — 7r/(2a) are pre- 
sented. 



located around T (see Fig. 7) shrink to this point when 
q x grows and disappear when q x becomes greater than 
k Q . 



The birefringence of the dispersion branches near the 
plasma frequency corresponds to the two signs in the 
right side of (38) and (39). The difference between equa- 
tions (38) and (39) is not significant if k « fco and 

In 4 the propagation of waves at the frequencies close to 
luq in triple wire medium has been numerically studied 
for the case when q x — (in- plane propagation). In 
that case the presence of a;— wires docs not influence the 
propagation characteristics, and our dispersion equation 
(37) is applicable. We compare the solution of (37) with 
results from 4 in order to validate our theory. 
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VI. QUASI-STATIC CASE 

Let us consider a double wire medium in a quasi-static 
case when |q|a <C 7r and ka <C n. Expanding trigonomet- 
ric functions in the dispersion equation for extraordinary 
waves (28) into Taylor series and keeping two first terms 
in those expansion, we receive to the following equation: 

(k 2 -q 2 y )(k 2 -q 2 z )[k 2 -k 2 (r y , a, c)-q 2 } [k 2 -k 2 (r z , a, b)-q 2 } 

= ll<£k%(r y) a,c)ko(r z ,a,b) (36) 

Let us consider a special case when r y = r z = r and 
a = b = c. In that case (36) could be simplified to the 
form: 

^jk 2 - ql^W^q 2 ^ 2 - k 2 -q 2 x -ql- q 2 z ) ± q v q z k 2 = 0, 

(37) 

where fco = ko(r, a, a). 

We can express q x from (37) in the form 

£ = k 2 - k 2 ± q ^ k \ - q 2 . (38) 

yjk 2 - q 2 ^k 2 - q 2 _ 

In 4 the following approximate dispersion equation was 
introduced under the conditions k ~ ko and qt = 
\J q 2 + q 2 "C ko (in our notations): 

q 2 x ^k 2 -kl±q y q z ~q 2 t . (39) 



FIG. 8: Dependence of the normalized wave number ka on tp 
angle near the "plasma" resonance of the wire medium with 
r = a/100, qt — O.ln/a (dashed line). Comparison with the 
exact result (solid line corresponds to the numerical data from 
Fig. 6, 4 ). 

In 4 one chose the following parameters: r = a/100 
(it corresponds to kga w 1.4), q t — 0.1w/a. One calcu- 
lates the medium dispersion for q x — at the frequencies 
k w fco- ln Fig- 6 of this work one shows the depen- 
dence of the normalized eigenfrequency ka on the angle 
ip. The angle (p is indicated in Fig. 4 (q y = q t sin tp and 
9z = qt cos ip). The plot ka vs. ip shown in Fig. 8 rep- 
resents the comparison of (37) with the numerical data 
from 4 . The upper dispersion branch corresponds to the 
case when before q y q z in (37) there is a plus sign, and the 
lower branch corresponds to the case when there is the 
minus sign. This plot illustrates the effect of the disper- 
sion branch birefringence near the "plasma frequency" 
of wire medium (see ellipses in Fig. 7). Equation (37) 
shows that this effect keeps also for q x ^ 0. Fig. 8 veri- 
fies that the quasi-static equation (37) is correct even at 
rather high frequencies (slightly higher than luq) outside 
the initial approximation ka -C n. 

Now, let us turn to the consideration of the effective 
permittivity of 2d WM. In the dyadic form the tensor of 
effective relative permittivity of arbitrary anisotropic di- 
electric media can be written as s = £ xx X- x. + s yy yo~yo + 
e zz zoZo. The dispersion equation of anisotropic dielectric 
(following from Maxwell's equation and from the defini- 
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tion of relative permittivity in terms of D = SqS ■ E) cor- 
responds to the zero- value determinant of the following 
system (f k 2 + qq q 2 I)~E = 0: 

(e xx k 2 -q%- q 2 z )E x + q x q y E y + q x q z E z = 
q x q y E x + (e yy k 2 -q 2 x - qf)E y + q y q z E z = (40) 
q x qzE x + q y q z E y + {e zz k 2 - q\ - q 2 )E z = 

This system helps to find the polarization of cigcnmodes 
when e is known. The study of eigenmodes polarization 
is going to be considered in the next paper. 

The dispersion equation has the following form: 

(q 2 y +q 2 z~ k 2 e xx ){q 2 x + q z - k 2 e yy ){q 2 x + q 2 y - k 2 e zz )- 



(q 



k 



\ 2 2 

c)q y q z 



{ql 



k 2 t 



2„2 



mi 



)q x q 



(q x +q y ~ k £zz)q x q y - 2q x q y q z = 0. 



(41) 



In 4 the following expressions have been heuristically 
introduced for components of the permittivity of 3d WM: 

ko(r x ,b,c) 



= 1 - 



k 2 



k"o{ r yi a i c ) 
k 2 -q^ 



= 1 



ko(r z ,a,b) 
k 2 -q 2 



(42) 



The effects of the low-frequency spatial dispersion are 
deemed to be described by terms q x , y , z in the denomina- 
tors of the components of e (see also in 3 ). 

It has been noticed in 4 that the expressions (42) and 
the dispersion equation (39) fit perfectly with the results 
of numerical simulations for ui s=a uj - Following (42), the 
components of e for triple WM are the permittivities of 
the three orthogonal simple wire media stretched along 
Cartesian axes. We have assumed that the same rule 
holds for 2d WM. In case of 2d wire medium there are 
no x— directed wires and e xx = 1. We have analytically 
verified that (36) exactly coincides with (41) if the effec- 
tive permittivity of a double WM takes the form: 



^double = X X + £yyVoyO + EzzZOZO, 



(43) 



where e yy and s 



are given by the relations (42). It 
should be noticed that the formula (43) has been ob- 
tained (very recently) for 2d WM by other authors as 
a result of a very complicated analytical-numerical ap- 
proach. ^From (42) it follows, that at every point of the 
central isofrcquency contour in Fig. 6 (where k < ko and 
qy,z < k) both components of the permittivity tensor e yy 
and e zz are negative. The propagation of such a wave 
(1^0 for it and the electric field can contain y— and 
z— components) is the spatial dispersion effect. 

We have also proved that the whole system (42) holds 
in our model of a triple wire media. The quasi-static 
analogue of (35) has the form 

(k 2 - q 2 x )(k 2 - q 2 )(k 2 - q 2 z )[k 2 - k 2 (r x ,b, c) - q 2 } X 



[k - k (r y ,a,c) - q }[k - k (r z ,a,b) - q } 



(k -q x )[k -k (r x ,b,c)-q }q y q z k (r v ,a,c)k (r z ,a,b) 



-(k -q y )[k -k (r y ,a,c)-q]q x q z k Q (r x ,b,c)k {r z ,a,b) 



(k -q z )[k -k (r z ,a,b)-q ]q x q y k (r x ,b,c)k (r y ,a,c) 



+^q x q y q z kQ{r x ,b,c)k Q {r y ,a,c)k Q {r Zl a,b) = 0. (44) 

It coincides with the dispersion equation (41) if the per- 
mittivity takes the form (see relations (42)): 



£triplo — SiiXoXq + £j,j,yoyo + ^zzZoZq. 



(45) 



We have thus verified, that dielectric permittivities for 
2d and 3d wire media in the form (43) and (45) sug- 
gested in 4 and 9 fit successfully in our theory. We have 
analytically verified, that the quasi-static analogues of 
dispersion equations (23) and (35) in the form (37) and 
(44) coincide with the dispersion equations of anisotropic 
dielectrics (43) and (45). 



VII. CONCLUSION 

In the present paper we have generalized a recently 
developed analytical theory of a simple wire media to the 
case of double wire media and obtained some results for 
triple WM. We have validated our theory by comparison 
with 4 and proved that the effective permittivity of 2d and 
3d WM introduced in 4 and 9 fits our dispersion equations 
fairly well. 

We have theoretically revealed the effects of low- 
frequency spatial dispersion for 2d WM such as: 

• Propagation of z— polarized TLM along y— wires 
is not suppressed by the presence of z— wires (the 
same is correct for the y— polarized TLM propagat- 
ing along z). 

• There are TLM which can exist in both y— and 
z— arrays simultaneously. These modes do not 
transport energy, since the directions of the cur- 
rents in wires are alternating along the a;— axis. 

• There are two propagating modes at low frequen- 
cies ui < ojq which are not TLM and not ordinary 
waves. One mode has non-zero electric field com- 
ponent in the plane (y — z) whereas both y— and 
z— components of the permittivity tensor are nega- 
tive. For the other one the isofrcquency contour is 
nearly hyperbolic. 

• Near the plasma frequency the two other waves ap- 
pear with crossing isofrcquency contours. 
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The materials under consideration could find various ap- 
plications due to the properties discussed in this paper. 
We would like to note especially such applications as cre- 
ation of low frequency super-prism and design of materi- 
als with negative refraction. Those properties of double 
WM will be discussed in the next paper. 

Finally, let us discuss the problem of the homogeniza- 
tion of WM. The equation (41) relates three unknown 
components of e, three components of the wave vector q 
and the frequency (or wave number k). The components 
of q are related through dispersion equation with k. It is 
clear that the problem of s has no unique solution in thus 
formulation. The same concerns 2d WM. Though (42) 
fits our dispersion equations this result is heuristic and 
the permittivity has been introduced and not derived. Is 
it reasonable to try to find other possible expressions for 
f? 

It is well-known that the effective material parameters 
of spatially dispersive media have another meaning that 
those of continuous media. The effective susceptibility of 
such media in presence of a point source depends on the 
source position and has nothing to do with the effective 



medium susceptibility for plane waves. The usual bound- 
ary conditions are not valid on the medium interface. So, 
the material parameters are not very helpful for solving 
the boundary problem for media with spatial dispersion. 
The goal of the homogenization of WM is modest: to 
describe the low-frequency propagating properties of an 
infinite medium in terms of those parameters. There- 
fore all we need is to introduce the permittivity which 
would 1) describe all effects we can reveal solving the 
correct (quasi-static) dispersion equation and 2) allow to 
find the polarization of eigenmodes correctly. For both 
2d and 3d WM the permittivity (42) comprises all the 
dispersion properties at low frequencies (ka < 1). As to 
eigenwaves polarization, the result (42) requests further 
studies which will be also presented in the next paper. 
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